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' Abstract 

O I We consider the stochastic background of gravitational waves produced by an early genera- 
^ '• 

^ ■ tion of Population III stars coupled with a normal mode of star formation at lower redshift. 

J> ! The computation is performed in the framework of hierarchical structure formation and is 

^ ' based on cosmic star formation histories constrained to reproduce the observed star for- 

b : 

- - . mation rate at redshift z < 6, the observed chemical abundances in damped Lyman alpha 
absorbers and in the intergalactic medium, and to allow for an early reionization of the Uni- 
verse at 2 ~ 10 — 20 as indicated by the first year results released by WMAP. We find that 
the normal mode of star formation produces a gravitational wave background which peaks 
at 300-500 Hz and is within LIGO III sensitivity. The Population III component peaks at 
lower frequencies (30-100 Hz depending on the model), and could be detected by LIGO HI 
as well as the planned BBO and DECIGO interferometers. 



1 Introduction 



In the last few years, there has been significant progress in our understanding of the early 
cosmic star formation history. The first year data obtained by WMAP indicates a large 
optical depth, implying that the universe became reionized at high redshift in the range 
11 < 2; < 30 at 95% CL IJ. To account for a period of early reionization, it has been argued 
that a generation of very massive stars preceded the oldest observed generation of Population 
II stars 12] • In addition, cosmic star formation rate (SFR) at z < 6, has been observed at 
levels significantly larger than the current rate [S]. Taken together, the evidence suggests 
that the distribution of the first stars (Population III) are described by a top-heavy initial 
mass function (IMF), formed in primordial metal-free structures with masses of order 10'' 
Mq. As these stars produced heavier elements, the universe achieved a critical metallicity 
(~ lO"'^ times solar metallicity) [T, at which point the massive mode of star formation 
yielded to a more normal distribution of stellar masses with a SFR peaked at 2; ~ 3. A more 
detailed understanding of the first epoch of star formation will rely on the phenomenological 
consequences of the models such as element enrichment and supernova rates El El • 

One consequence of this new view of star formation is an enhanced rate of core collapse 
supernovae. The resultant relic neutrino background was investigated in Ref. [Hlini- In each 
core collapse supernova explosion, the bulk of the energy released is in the form of neutrinos 
which, because they are weakly interacting, retain information about their origins. Although 
the neutrino background from a massive mode of star formation (Pop III) at early times is 
not likely to be detected due to its redshifted spectrum, the prospects for observation of the 
spectrum produced by the normal mode of star formation (Pop II) in the near future are 
good. 

One possible probe of the massive mode, however, is the stochastic background of gravi- 
tational waves produced by cosmological core collapse supernovae, which we consider here. 
Although only a small fraction of the total energy of core collapse is emitted in gravitational 
radiation, the improved sensitivities by way of correlation of currently operating ground- 
based interferometers GEO600 [IHI, LIGO II and III [11, TAMA P], and VIRGO 
and of future space-based antennas BBO ^5^, LISA [15j, DECIGO make the positive 
detection of the accumulated gravitational wave background plausible. 

The gravitational wave background from core collapse supernovae resulting in black holes 
has been calculated in Refs. ^Tj , with estimates of the peak of the differential energy 
density spectrum of flcwh'^ = 10~^^ — few x 10~^ reaching its maximum value at frequencies 
anywhere from a few hundred Hertz to a few thousand Hertz. The calculation has also been 
made specifically for Population HI supernovae resulting in black holes in Ref. They 
estimate the spectrum to peak at flcwh"^ ~ 10~^ at a frequency of O(IOO) Hz. More recently, 
the spectra from both a normal mode of star formation, in which all stars collapse to form 
neutron stars, and a Population HI mode, in which all stars collapse to form black holes, was 
calculated in Ref. [201 • They find both peaks to be located at roughly Qc^yh'^ ^ few x 10^^^ 
in the most optimistic case, with the peak frequency dependent on the redshift range over 
which gravitational collapse occurs. As these spectra are highly model dependent, and given 
that previous estimates of the differential closure density span four orders of magnitude in the 
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range where detection may soon be possible, it is now important to examine the sensitivity 
of a detectable signal to the star formation history. 

Here, we incorporate fully developed chemical evolution models which trace the history of 
pre-galactic structures as well as the IGM and are based on a ACDM cosmology with a Press- 
Schechter model of hierarchical structure formation We adopt the chemical evolution 
models of Daigne et al. and consider several bimodal star formation histories, each with a 
normal component of star formation as well as a massive component describing Population 
III stars. Given an IMF and a respective SFR, we calculate the expected gravitational wave 
background and compare this result with detector sensitivities. 



2 Calculation of the Gravitational Wave Background 

Gravitational waves can be characterized by a frequency, /, and an amplitude, /i, which 
is defined by the degree of quadrupole anisotropy and strength of the source. Using an 
amplitude determined from the simulation of a 15 M0 star from Ref. [22], Buonanno et 
al. [21? made a generalization to larger stars by using the function's dependence on the 
anisotropy and neutrino luminosity during collapse. The shape of the gravitational wave 
spectrum can then be described by the dimensionless quantity 
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where h is the Fourier transform of h, Gn is Newton's constant, Ei, is the total energy 
emitted in neutrinos, (q) is the average value of the anisotropy parameter, q, defined in 
Ref. [23], and D is the distance to a typical supernova. Although it is necessary to know the 
distance in order to quantify the amplitude, we will see that the accumulated energy density 
does not depend on this parameter. Once the structure of the amplitude is imposed, the 
constants a and b determine the specific spectral shape. By roughly reproducing the source 
spectrum from simulation model sl5r in Ref. one obtains a ^ 200 Hz. and b ^ 300 Hz. 
The sensitivity to these choices will be discussed in section 

A stochastic background of gravitational waves with energy density pcw and frequency 
/ is best described by the differential closure density parameter |2_4|, 
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with Pc = 3Hq/8ttGn the critical density. Given a star formation history, consisting of a 
star formation rate per comoving volume (SFR), ip{t), and an initial mass function (IMF), 
4>{m), and the gravitational wave amplitude, the stochastic background of gravitational 
waves produced by cosmological core collapse supernovae is 
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where /' is the frequency at emission, related to the observed frequency, /, by /' = /(I + z), 
Zi is the initial redshift at which supernovae begin to occur, Mmin and Mmax are the minimum 
and maximum masses in each model for which supernovae occur, and r(m) is the lifetime of 
a star of mass m ^2iij. Using Eq. ^ one obtains 
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The parameter (g) must be determined from simulations. As in Ref. [201; istke (q) = 
0.45% for progenitors that collapse to neutron stars, as determined by simulation in [22j. 
The total energy emitted in gravitational waves can be described in terms of the amplitude 

as 
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The efficiency of gravitational wave production, e, is defined in terms of the remnant mass, 
Mr, by 

Egw = eMrC^. (6) 

Using = 3 X 10^^ ergs, we determine the efficiency to be e = 1.5 x 10^'' when a neutron 
star remnant is produced. 

When collapse proceeds to form a black hole, we assume, as in Ref. jSHj, an efficiency 
of gravitational wave production found by Fryer et a/, of e = 2 x 10"^ [25J. One can then 
determine the quantity {q)Ei, for a star of mass m by assuming only a spectral shape (here, 
Eq. and an efficiency of gravitational wave production. In this case it is necessary to 
specify the mass of the remnant, but not the neutrino luminosity. Note that the efficiency 
of gravitational wave production is much greater when a black hole rather than a neutron 
star is produced as a remnant. 

One interesting consequence of choosing the efficiency of gravitational wave production 
to be constant for stars that collapse to black holes is that the anisotropy parameter actually 
decreases for larger stars, as can be seen in Fig.HJ This is because the quantity (g)-E'jy oc VM^, 
but the energy emitted in neutrinos itself, E^, oc Mr, resulting in (g) oc l/^/M^. The values 
obtained for (q) are not unreasonable for larger stars in Models 1 and 2b, however it has 
been noted that one might expect a larger degree of anisotropy due to rotation and/or 
violent explosions [20] • Very little is known about the anisotropy, however there is evidence 
that stars that collapse to black holes should be more efficient emitters of gravitational 
waves. Stark and Piran found a maximum efficiency of e < 7 x 10^"^ for an axisymmetric 
collapse resulting in a black hole [2^] , and more recently. Fryer et al. obtained the efficiency 
of 2 X 10~^ for a 100 M© black hole remnant, which we adopt. Note that since ^gw oc e, 
if the efficiency is actually closer to Stark &; Piran's maximum value, the observed energy 
density in gravitational waves may be as much as 35 times larger than that shown in our 
results below, and even more if the collapse is less symmetric. While it is possible to fix the 
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Figure 1: Anisotropy parameter, (g), as a function of the mass of the progenitor. The 
energy emitted in neutrinos is taken to be ii^j^ = 3 x 10^^ ergs for 8 Mq < m < 30 Mq 
and = ^M^c^ for stars with m > 30 Mq that collapse to black holes. For stars with 
30 Mq < m < 140 Mq, the mass of the remnant is taken to be the mass of the helium 
core before collapse, Mhs = 13/24 ■ (m — 20Mq) [27]. In the case of the larger stars with 
m > 270 Mq, Mr = m. Since (q) oc l/E^, this is a minimum average value for the anisotropy 
parameter for these larger stars. 



anisotropy parameter to be a function of mass, there is no data to support any particular 
functional dependence. Note that the gap in Fig. ^ falls in the range where pair-instability 
supernovae are important. The gap reflects a lack of knowledge concerning the production 
of gravitational waves for these stars. 

In all of our computations, we assume a flat ACDM cosmology with 
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where = 0.73, Vt^ = 0.27, and h = 0.71 p. 



3 Star Formation Models 

The cosmic star formation histories considered here have been adopted from the detailed 
model of chemical evolution in Ref. [7j. The models are described by a bimodal birthrate 
function of the form 

B{m,t,Z) = Mm)Mt) + <P2MMZ) (8) 

where 0i(2) is the IMF of the normal (massive) component of star formation and ipi{2) is the 
respective SFR. Z is the metallicity. The normal mode contains stars with mass between 0.1 
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Mq and 100 Mq and has a SFR which peaks at 2; ~ 3. The massive component dominates 
at high redshift. The IMF of both modes is taken to be a power law with a near Salpeter 
slope so that, 

0i(m) oc m-(^+^) (9) 
with X = 1.3. Each IMF is normalized independently by 



dm m(j)i{m) = 1 , (10) 



differing only in the specific mass range of each model. Here we consider two different mass 
ranges for the massive mode, although three mass ranges are presented in Ref. 0, as will 
be discussed presently. Both the normal and massive components can contribute to the 
chemical enrichment of galaxy forming structures and the IGM, though the normal mode is 
not sufficient for accounting for the early reionization of the IGM 

Here, we restrict our attention to the best fit hierarchical model in [Zj in which the 
minimum mass for star formation is 10^ Mq. The normal mode SFR is given by 

^i(t) = UiMstructexp (-t/ri) , (11) 

where ti = 2.8 Gyr is a characteristic timescale and ui = 0.2 Gyr~^ governs the efficiency of 
the star formation. In contrast, the massive mode SFR is defined by 

tp2{t) = I^2MisMexp (-ZiGM/^crit) , (12) 

with z/2 maximized to achieve early reionization without the overproduction of metals or the 
over- consumption of gas. We adopt Z^ntl^Q = 10""^ . 

We consider three different models, labeled Models 1, 2a, and 2b to describe the massive 
mode. They are distinguished by their respective stellar mass ranges. In Model 1, the IMF 
is defined for stars with masses, 40 Mq < m < 100 Mq. All of these stars die in core collapse 
supernovae leaving a black hole remnant. Model 2a is described by very massive stars which 
become pair instability supernovae. The IMF is defined for 140 Mq < m < 260 Mq. Finally, 
the most massive stars are considered in Model 2b and fall in the range 270 Mq < m < 500 
Mq, with the SFR as in Model 1. These stars entirely collapse into black holes and do 
not contribute to the chemical enrichment of either the ISM or IGM. The coefficient of star 
formation, 1^2 is 80, 40, and 10 Gyr~^ for Models 1, 2a, and 2b respectively. In each case, 
star formation begins at very high redshift {z ~ 30) but peaks at redshifts between 10 and 
15, depending on the model. Note that the absolute value of the SFR depends not only on 
1^2, but also on the efficiency of outflow. See [Zj for details. 

In Figure 12 we show the SFR, for Models 1, 2a, and 2b (including the normal mode). 
Also shown by the dotted curve is the SFR for the massive mode alone in Model 1. For 
comparative purposes, we also consider an example of a model (with an IMF as in Model 
1) in which the massive mode occurs as a rapid burst at 2; = 16 designated as Model le. 
This is also shown in Figure |21 compared with the analogous result for Model 1. As one can 
see, the massive burst SFR is significantly larger than the model considered above for the 
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Figure 2: The star formation rate for Models 1, 2a, and 2b as labeled by their respective 
mass ranges, and a rapid burst model. The dashed line in the upper panel shows the SFR 
of the massive mode of Model 1. 

short duration of the burst, while at lower redshifts the SFR, which is determined by vi is 
nearly identical. Because of our lack of understanding of gravitational wave production in 
pair-instability supernovae, we will not consider Model 2a any further. 

The rate of core collapse supernovae can be calculated directly in terms of the IMF and 



Jmax(8Af0,m,„i„(t)) 

where mmin{t) is the minimum mass of a star with lifetime less than t. The differential energy 
density parameter for each model is calculated using Eq. HI 

4 Results 

The differential energy density in gravitational waves from core collapse supernovae in Model 
is shown in Fig. El These stars collapse to form either neutron stars or black holes with a 
mass equal to the mass of the progenitor's helium core. The spectrum peaks at a frequency 
of about 360 Hz, with VLcwh^ = 3.5 x 10~^°. This spectrum should be similar to the upper 
curve in Fig. 4 of Ref. which peaks at roughly 300 Hz but with a maximum value 
about two orders of magnitude lower than that in Fig. El One can see from Fig. |31 that our 
supernova rate lies below theirs at high redshift and they are comparable at later times. 
A large supernova rate at high redshift serves to broaden the gravitational wave spectrum, 
but in this case it is not the cause of the disparity in peak heights, which is attributable 
primarily to differences in the IMFs. Their spectra were calculated under the assumption 
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Figure 3: Gravitational Wave background from core collapse supernovae in Model 0. 

that all progenitors collapse to form neutron stars, whereas our Model also includes stars 
that collapse to black holes. As can be seen from Eq. the number of stars collapsing to 
black holes in Model is suppressed by an IMF that favors low masses, but stars emit ~ 3 
orders of magnitude more energy in gravitational waves when the collapse yields a black 
hole than when it yields a neutron star due to the larger efficiency and remnant mass. As a 
result, the gravitational wave background for the normal mode of star formation is enhanced 
by two orders of magnitude over that in Buonanno et al. pOj. 

The background of gravitational waves from stars that collapse to black holes was cal- 
culated by Araujo et al. [TH] assuming a Springel & Hernquist j2H] model of star forma- 
tion. They find that for a Salpeter IMF defined for stars up to 125 M© and all stars 
with 25 < m < 125 collapsing to black holes with Mr = |m, the spectrum peaks at 
^Gwh"^ ~ 5 X 10^^ at / ~ 200 Hz. This amounts to an expected background about an order 
of magnitude larger than that from Model 0. The Springel & Hernquist SNR is actually 
smaller than any normal mode SNR considered here by a factor of ~ 8 as one can see in 
Fig. 0J but they assume a much larger efficiency found in Ref. [^Hl of e = 7 x 10~^ and this 
remnant mass fraction is actually larger than any found in our Model 0. 

Figure IHl shows the gravitational wave background from the massive modes of Models 1, 
le, and 2b. In each case, the SFR has a different shape and amplitude. One can see that 
the spectrum from the massive mode of Model 2b is quite similar to the sharply peaked SFR 
case of Model le, with the peak location shifted due to star formation reaching maximum 
at different redshifts. As a burst, star formation in Model le occurs in a very limited 
redshift range around z = 16. As a consequence, the frequency of gravitational waves 
does not extend much past 100 Hz. In contrast, star formation in Model 2b does extend 
down to lower redshifts and hence we expect the spectrum to extend to higher frequencies. 
Similarly, the gravitational wave spectrum from the massive mode of Model 1 extends to 
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redshift 

Figure 4: Rates of core collapse supernovae considered in the calculation of the gravitational 
wave background for Model (solid blue), from Buonanno et al. (dotted red) ^20j, and using 
the Springel and Hernquist SFR from Araujo et al. (dashed black) [THj . 

still higher frequencies due to significant star formation at redshifts between 3 and 5. The 
fact that maximum values of the energy density parameter span two orders of magnitude is 
attributable primarily to the differing SFRs seen in Fig. but also to the differing remnant 
masses in Model 1 and 2b. Stars in Model 1 have much smaller initial masses and collapse to 
form black holes with a mass equal to the progenitor's helium core, whereas stars in Model 
2b, which are larger to begin with, collapse entirely to black holes. This entails larger bulk 
motions of matter, with VLqw oc M^, as can be seen from Eqs. 

The spectrum from supernovae in Model 2b, which includes 300 Mq progenitors, could 
be compared with the curves in Figure 8 of Ref. [201 ■ Although the frequencies at which the 
peaks occur are similar, the energy density obtained in our calculation is larger. The differing 
peak heights are attributed predominantly to the difference in the assumed supernova rates 
for the Pop III modes. In Model 2b, the baryon fraction in Pop III stars is //// = 7.0 x 10^^. 
This is much larger than the maximum value of //// = 10"'^ assumed in [20]; yielding VLcwh^ 
approaching 10~^^. The supernova rate, and therefore the energy density in gravitational 
waves, scales linearly with ////. In addition, we assume total collapse of massive stars 
leaving a remnant mass equal to the progenitor mass whereas Buonanno et al. assume that 
300 Mq stars in collapse to form 100 M© black holes, and VLcwh'^ is linear in the remnant 
mass. Finally, Model 2b, is not a delta function of 300 M0 stars, but rather a power law 
distribution containing stars up to 500 M©. These effects combine to enhance our peak 
height by a factor of about 300. It is more difficult to compare Models 1 and le, as the IMF 
contains stars with masses less than 100 Mq. Our fraction of baryonic matter in Population 
III stars is different for each model; 3.3 x 10^'^ for Model 1 and 1.2 x 10^^ for Model le. 

The total gravitational wave background for Models 1, le, and 2b, including the normal 
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Figure 5: Gravitational wave background from Population III core collapse supernovae in 
Models 1, le, and 2b. 
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Figure 6: Gravitational wave background from core collapse supernovae in Models 1, le, and 
2b for both the normal and massive mode combined. 

mode, is shown in Figure IHl In the spectrum for Models le and 2b, one can clearly see the 
peak from the normal mode near 360 Hz and that from the massive mode at lower frequency. 
The massive mode contribution to the total background spectrum in Model 1 is not large 
enough to be seen here. Consequently, the spectra for Model alone and Model 1 (including 
the normal mode) are identical. However, as discussed in Sectional Model alone could not 
provide a sufficient flux of ionizing photons for reionization at high redshift. 

As was mentioned in Section the shape of the spectrum of gravitational waves depends 
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on the parameters a and h in Equation IHl While any variation in these parameters can have 
a dramatic effect on the resultant spectrum when the anisotropy parameter and neutrino 
luminosity are specified, the effects are strongly suppressed by our normalization scheme in 
which we specify the efficiency of gravitational wave production, e, for stars that collapse to 
black holes. The effect of varying the parameters a and h in the spectrum for Model 2b is 
shown in Figure H 




0.1 1 10 100 1000 



f(Hz) 

Figure 7: Total gravitational wave background from core collapse supernovae in Model 2b 
including both the normal and massive modes. The three curves correspond to different 
values for the parameters a and h. Our nominal case of a = 200 Hz and h = 300 Hz is shown 
as the solid curve. The dashed (dotted) curve corresponds to a = b = 250 Hz (a = 150 Hz 
and b = 350 Hz). 



5 Detection 

As the backgrounds calculated here are the collective result of the collapse of all supernovae 
in the models considered, it is important to know whether the background is continuous. 
The ratio of the duration of each burst to the time between successive bursts is called the 
duty cycle, given by 

DC= r dRsN^{l + z) (14) 
Jq 

where dRsN is the differential rate of core collapse supernovae as seen from Earth and Ar 
is the average time duration of a single burst of gravitational wave emission [22] • One can 
approximate the duty cycle by choosing a fixed value for Ar or by noting that the time 
duration of the burst corresponds to a frequency at emission, so Ar(l + z) ~ 1//, where / 
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is the observed frequency. Then the duty cycle can also be approximated by 



BC . i /" Ms. ^ i r- (15) 

f Jo f Jo Jm^,„ 1 + ^ dz 

where dV/dz = 47Tr'^{z)/H is the comoving volume element and dr = [1 + z)dt The 
duty cycles approximated using these two methods are given in Table In column 2, the 
duration of the burst was assumed to be roughly the same for all supernovae, At ~ 1 ms. 
As DC < 1, the background is not expected to be continuous. Using the method described 
by Eq. E[ we find that the background is expected to be continuous for low frequencies 
(/ ~ 10 Hz.), but that in each case individual bursts will be distinct (shot noise) throughout 
the peak of the spectrum. 



Model 


DC(Ar = 1 ms) 


-DC x/ 





0.18 


37 


1 


0.23 


41 


le 


0.47 


54 


2b 


0.34 


61 



Table 1: The duty cycles calculated for each of the models considered here. All include both 
the normal and massive modes, except Model 0, which is the normal mode only. 



Finally we comment on the potential for current and future detectors to observe the 
gravitational signal from the first stars. Figure |H1 shows the sensitivities of the LIGO II 
and III correlated ground based detectors and proposed sensitivities for LISA, the Big Bang 
Observer (BBO) and the Decihertz Interferometer Gravitational Wave Observatory (DE- 
CIGO). LISA, intended to search for gravitational waves from merging binaries as well as 
any primordial background, is designed to be sensitive to low frequencies < 0.1 Hz. Un- 
fortunately, comparison of LISA's sensitivity curve with the dashed lines representing the 
expected backgrounds for Models le and 2b show that LISA will be incapable of detecting 
the backgrounds predicted here. DECIGO and, as a follow on to LISA, BBO have been 
proposed to search for the inflationary gravitational wave background at higher frequencies. 
Both antennas will have sufficient sensitivity to detect the gravitational wave background 
from cosmological supernovae at frequencies / ^ 10 Hz. Based on CMB temperature fluc- 
tuations, the upper limit to the gravitational wave background from inflation was estimated 
to be Qcwh'^ = 10~^^ which may be completely masked by that from stellar collapse 
down to frequencies of 10~^-10~^ Hz. 

With a lead time of more than 20 years for BBO (and similarly for DECIGO), LIGO 
correlated detectors are the most immediate promising prospect for detection. Although 
LIGO II will not have the requisite sensitivity to detect the background from any model 
considered here, LIGO III correlated detectors will be able to probe the frequency range 
/ ~ 10^ - few xlO^ Hz. to sufficiently small values of the energy density that the peak from 
the massive mode may be observed. When considering the frequency ranges covered by both 
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LIGO III and BBO/DECIGO, it is possible that the entire spectrum below a few hundred 
Hz. may be observed in the near future. 
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Figure 8: Sensitivities of proposed and operating gravitational wave detectors BBO, DE- 
CIGO, LIGO II and III (correlated), and LISA expressed in terms of energy density 
The dashed lines are the expected gravitational wave backgrounds from Models le and 2b 
as in Figure ini 



6 Conclusions 

We have calculated the gravitational wave background from three different star formation 
histories that reproduce the observed chemical abundances and reionize the universe at high 
redshift. Each star formation history consists of a Population III mode of star formation 
coupled to a normal mode that reproduces the observed SFR for z < 6. The gravitational 
wave background was calculated assuming the amplitude given in Eq. ^ We obtained three 
distinct gravitational wave background spectra. In each case, there is a peak due to the 
normal mode near / = 360 Hz. as well as a peak due to the massive mode at lower frequency. 
The location of the massive mode peak depends on the redshift range over which Population 
III star formation occurred; the earlier star formation reached maximum, the smaller the 
peak frequency. 

The gravitational wave background from core collapse supernovae in our models is found 
to be large enough that it should be detected by the next generation of space-based laser 
interferometers, BBO and DECIGO, if not sooner by LIGO. The background will constitute 
a shot noise signal in the frequency range accessible with LIGO, however we expect the signal 
to be continuous for frequencies that will be probed by BBO and DECIGO. Within the next 
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few decades, a detection of the gravitational wave background from cosmological supernovae 
will provide valuable information about the history of structure formation in the universe. 
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